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Abstract. We consider the possibility to extract spins that are generated by an 
electric current in a two-dimensional electron gas with Rashba-Dresselhaus spin- 
orbit interaction (R2DEG) in the Hall geometry. To this end, we discuss boundary 
conditions for the spin accumulations between a spin-orbit coupled region and contact 
without spin-orbit coupling, i.e. a normal two-dimensional electron gas (2DEG). We 
demonstrate that in contrast to contacts that extend along the whole sample, a spin 
^ . accumulation can diffuse into the normal region through finite contacts and detected by 

£j \ e.g. ferromagnets. For an impedance-matched narrow contact the spin accumulation 

in the 2DEG is equal to the current induced spin accumulation in the bulk of R2DEG 
up to a geometry-dependent numerical factor. 
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1. Introduction 

In recent years, there has been an increasing impetus towards generating and detecting 
spin accumulations and spin currents in nonmagnetic systems. Conventional means 
of achieving this goal are to use ferromagnets and magnetic fields to inject and/or 
detect spins pQ. Recently, spin generation based on two related effects, current induced 
spin accumulation [21 El S] and current induced transverse spin current [6] (known 
as the spin Hall effect), has attracted considerable attention. In Ref. [6], the spin 
Hall effect was caused by the spin-orbit (SO) interaction of impurities and the effect 
is then called "extrinsic". The "intrinsic" SHE caused by a band structure with 
SO-induced spin splittings was proposed by Sinova et al. [9J for the R2DEG and 
Murakami et al. [10] for the hole gas in bulk III-V semiconductors with significant 
SO interaction. After an initial controversy, it is now generally agreed that in the 
diffuse regime the SHE vanishes in the bulk of a 2DEG with k- linear (Rashba and/or 
Dresselhaus) SO coupling [HI [TJl [13], but remains finite for extrinsic SO coupling, 
intrinsic SO coupling in two-dimensional hole systems, and near the edges of a finite 
diffusive R2DEG [T2l [14]. The spin Hall effect has been observed in semiconductor 
electron [16] and hole [TF] systems by the detection of edge spin accumulations 
with optical methods, and in metals by the electrical detection of spin currents via 
ferromagnetic leads [IE]. Although initial theoretical investigations of the SHE and 
current-induced spin accumulation has been on bulk disordered conductors using Kubo, 
Keldysh or Boltzmann formalism [7] El HQl HH H21 HSl HH HHl [20] , it is now understood 
that the bulk conductivity is not necessarily related to experimentally relevant quantities 
such as local spin accumulations probed by local optical or electrical probes. In this 
respect, a more local approach based on spin diffusion equations is advantageous [T2| [13] . 
However, spin diffusion equations have to be supplemented by suitable boundary 
conditions that have observable consequences. There have been many proposals in 
that direction [TJl [22], [211 1231 [MJ (251 [26], but a consensus has not been reached so far. 

Here, we focus on the boundary conditions between a (half infinite) 2DEG with 
finite Rashba type spin orbit coupling (R2DEG) and a (half infinite) 2DEG without 
spin-orbit coupling connected by a contact that is narrow on the scale of the system, but 
wider than the mean free path. Such a boundary has been considered by Refs. [23, 26J, 
but for an infintely wide contact region, for which it could be shown that no spin 
accumulation could diffuse into the 2DEG [26] . We shall show below, however, that for 
a narrow (as opposite to wide) contact, the spin accumulation in the 2DEG is equal to 
the bulk value of the spin accumulation in R2DEG up to a numerical constant which 
depends on the geometry that is smaller but can be of the order of unity. These results 
prove that current induced spins can be extracted to a region with small spin-orbit 
coupling in which the spin lifetime is very long and used for spintronics applications, 
thus confirming our previous results [T4"] . 

This article is organized as follows: we define our model and derive spin diffusion 
equations in section [2j In section [3J we first recapitulate the symmetry relations for 
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conductances with respect to measuring the spin accumulation in a normal region 
with ferromagnetic leads. Next we apply these relations to demonstrate that the spin 
accumulation from the R2DEG can be extracted into a 2DEG region. In section HI we 
focus on a model for a small contact between the R2DEG and the 2DEG and solve it to 
demonstrate the principle of spin extraction to a region with vanishing SO interaction. 
The numerical simulations for the diffuse R2DEG — 2DEG heterostructure are reported 
in section [51 



2. Spin diffusion equations in a 2D electron gas with Rashba spin-orbit 
coupling 

In this paper we focus on a disordered finite size 2DEG with Rashba type spin-orbit 
coupling, noting that the effects of a significant Dresselhaus term can be included 
straightforwardly. Throughout the paper we shall assume that all length scales of this 
finite region are much larger than the elastic mean free path such that spin transport 
is governed by diffusion equations [12J, [13] is valid. In this section, we proceed to derive 
these spin diffusion equations for later convenience. 

In 2 x 2 spin space, our system is defined by the Hamiltonian: 

P 2 

H = £- + ap-((TXz) + C/(x) + V(x) (1) 
2m 

where x and p are the (two-dimensional) position and momentum operators, 
respectively, a is the vector of Pauli spin matrices (the 2x2 unit vector is implied with 
scalars), z is the unit vector normal to the 2D plane, and a parameterizes the strength 
of the SO interaction that can be position dependent e.g. due to local external gates, 
and V(x) = z~2i=i <K X ~~ -^i) the impurity potential, modelled by N impurity centers 
located at points {Xi}, which for the sake of simplicity we assume to be spherically 
symmetric, U(x) is a smooth potential that confines the system to a finite region but 
allows a few openings to reservoirs. 

2.1. Rashba Green function 

Our starting point is the impurity averaged Green's function G(k) = (h 2 k 2 /2m + harj ■ 
k — E — ih/r)^ 1 , where r\ = z x er, and r is the momentum lifetime. In terms of its 
components, G{k) is given by 

2m _„ ,1/1 1 
- G ( k ) = o T^^T + 



h 2 w 2 V k 2 - kl k 2 - k 2 _ 



fcqk ■ T) k 2 j2 - ., ■ t , , 

" l " p _ p I p _ p p _ p J ' v z J 



where k\ = k 2 F + k 2 j2±k a ^k 2 F + k 2 jA + 2mi/{rh), k a = 2ma/Hand k F = ^2mE F /h 2 . 
The real space Green function is then obtained by a Fourier transform: 

G(x, E F ) = ^ - l -(H (k + x) + H (k_x)) - J^_(H (k + x) - H (k-x)) 
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(3) 



where x = |x|. We note that we only need the large kpx asymptotics of G(x), because 
we are interested in dilute disorder. The conventional approximation [27] is to expand 
G(x) to leading order in l/{kpr) and k a jkp\ 



G(x, E F ) « -^J- e ft^M2i e -«.^ ( (4) 

where Z = kkpr/m . This level of approximation is sufficient for most spin-orbit related 
applications such as the calculation of Dyakonov-Perel spin relaxation, spin precession, 
weak antilocalization etc. However in order to study current induced spin accumulation 
and SHE in diffusive systems, it is necessary to go to higher order in ma/hkp and 
1/(/cfx). With these correction terms the asymptotic Green function becomes: 



3i 



-e 



ifc a x-Tj/2 _| ^ik a x/2 _|_ e ~ik a x/2^ 



8k F x ' 8k F x v 

where x = x/x In the next subsection, we will use this expression to derive spin diffusion 
equations for a R2DEG. 

2.2. Diffusion equation 

We first focus on the equation of motion of the density matrix with coherent spin 
components. It can be shown that in the limit Epr/h ^> 1, the energy resolved density 
matrix satisfies the following equation [121 EES [23j 121] : 

p a (x,u) = — — / dV/C ac (x,x';^)p c (x',^), (6) 

27TUT J 

where p a = tr(pa a ), summation over repeated indices is implied, and v is the density of 
states and 

/C a ,(x,x» = ^Tr( K a a G R (x,x';E + u;)a b G A (x',x ] E)y (7) 

Multiplying p(E) with the density of states and integrating over energy, we obtain the 
densities and polarizations, whereas accumulations are obtained by directly integrating 
over energy. The diffusion equation is obtained by expanding Eq. [7J to second order in 
spatial gradients. In a homogeneously disordered system we have: 

p a (x) = — — / d 2 rfC ac {r)p c {Y + x) 
Zttut J 

~ -— / d 2 r/C ac (r)(p c (x) + r ■ Vp c (x) + r i r i 9 i 9 i p c (x), (8) 

LlXVT J 

where p a (x) = p a (x; 0). We now use the asymptotic expression Eq. (jSJ) for the Green's 
function and insert the resulting expression in to Eq. (IS]). The spatial integrals are 
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elementary and lead to the following equations for the vector components of the density 
matrix, Sj = pi/ 2 and n = p : 

DV 2 n - 4K S _ C (V x s) z =0 (9) 



DV 2 s 3 - 2K P (V • s) = — (10) 



DV 2 s + 2K p Vs 3 - K s . c (z x V)n = — (11) 

T s 

Here D = t s = r(l + 4£ 2 )/2£ 2 (the Dyakonov-Perel spin relaxation time), 

K s _ c = a!; 2 /(I + 4£ 2 ), K p = Hkp£/m(l + 4£ 2 ) 2 and £ = appr/h. A similar expansion 
for the spin current, this time to first order in the spatial gradients, produces the analog 
of Fick's law for spin diffusion: 

i) = r~T772 ( ^3 (sj - e jm3 —V m n) - 5 i:j s 3 ) - vDVjSi. (12) 



T 8 
S 



I _|_ 4^2 \ 10 \ j J ' no 2 

When supplied with suitable boundary conditions the diffusion equations (T9WTTT) and the 
spin current expression (TI2"]) can be solved to obtain all spin and charge conductances. 
Here, we are mainly interested in the boundary between a R2DEG and a 2DEG (for hard 
wall boundary conditions see Refs. [211 [231 El] ) . I n this case, the boundary conditions 
require the continuity of the spin current [HI [26] 

-^yn.^-^.fnxVln)-^) -vDn-Vsf. = uDn • Vsf ,(13) 
1 + 44 VV2 / / 

where s R and are the spin accumulations in the R2DEG and the 2DEG respectively, 
and n is the unit normal vector at the interface. A common choice for the matching 
condition for the spin accumulation at the interface is to assume that the spin 
accumulations are continuous (see e.g. Ref. [T]): 

L- < 14 > 

This condition has been criticized recently in Ref. [26] in which it was demonstrated 
that for an infinite interface with a constant electric field parallel to it: 

ar 



s N 



+ _ ^ _n ( n ■ (z x Vn) 





= 8^ 


)) 










(15) 



We first note that when the charge current is perpendicular to the interface, such 
as for a two-probe configuration [42], these two boundary conditions agree and no 
controversy exists. However, for an infinite interface where the charge current density 
is homogeneous, the difference between these two boundary conditions is drastic: if 
Eq. (1141) is valid, a current induced spin accumulation diffuses into the 2DEG. On the 
other hand, if Eq. f[T5l) is valid, the spin accumulation vanishes in the 2DEG. We solve 
this conundrum below by showing that for a contact smaller than the spin relaxation 
length (as assumed in Ref. [14"]). the two boundary conditions lead to results that agree 
up to a numerical factor of the order of unity. We therefore conclude that it is possible 
to extract spin accumulation to the 2DEG and detect it with a ferromagnet. 
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Figure 1. Setup for detection of current induced spins 



3. Onsager's relations and the spin boundary conditions 

In this section we provide a general symmetry argument based on Onsager's relations, 
that that proves viability of electric detection of the SHE and the current induced 
spin accumulation by finite size contacts. Let us start by addressing the symmetry 
properties of multiprobe conductances relevant for the combination of a spin-orbit 
coupled region with a ferromagnet via a normal region (Fig. [T]), using Onsager's 
relations [281 ISS [301 EU [32]. We are particularly interested in the setup shown in 
Fig. [TJ The configuration in Fig. []Ji is designed to measure the spin accumulation 
in the 2DEG injected from the neighbouring R2DEG. The voltage signal V directly 
observes boundary conditions between R2DEG and 2DEG when the charge current 
is parallel to the boundary. The setup in Fig. [lb, on the other hand, measures how 
much spin is injected into the R2DEG from the ferromagnet through the 2DEG. V 
measures directly the spin boundary conditions for a charge current perpendicular to 
the boundary. Onsager relations relate these two conductances, enabling us to relate 
the boundary conditions when the current is parallel or perpendicular to the boundary. 



3.1. Onsager's relations 

A generic SO-coupling operator consists of combinations of velocity and spin operators 
that are invariant under time reversal. When the spin-orbit coupled region is brought 
into contact with a Ferromagnetic region, the Hamiltonian of the combined system 
has the symmetry TH{rn)T~ l = H(—m), where m is a unit vector in the direction 
of the magnetization of the ferromagnet and T the time-reversal operator. We now 
focus on a the specific four-probe setups in Fig. [2] (for a more general discussion see 
Ref. [32]). The currents in the leads and the respective chemical potentials of the 
reservoirs are related in linear response as J, = YlijGijUj- We now use the Landauer- 
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Biittiker formalism to obtain Gij. The scattering matrix for the spin orbit (SO) 
coupled region and the ferromagnetic region is given respectively by Sso an d S m . The 
symmetry properties of these matrices are self-duality (reflecting the presence of spin- 
orbit coupling) Sso — S2>S'J S 2 , and S m = £ 2 ST m £2, where S 2 is block diagonal in the 
Pauli matrix a y [33]. We are interested in the block structure of Sso singling out lead 
3 combining the SO and F regions: 




Sso = I r T I (16) 
so 

where the matrix rso includes all reflections and transmissions that begin and end in 

the leads 1, 2 and 4. Using the rules for combining S-matrices, we obtain the joint 
S- matrix of the combined SO|F region: 

t = t m [l - r' so r m )~H so (17) 

t' = t[ so [l - r m r' so }-H' m (18) 

r = r so + t'so r m[l - r' so r m )~H S o (19) 

soC (20) 




' / ii r-i / „ i— 1 / 



r 



S = , (21) 



Using these rules we obtain the symmetries of the combined S matrix: E 2 t T (m)X 2 = 
t'{— m) and S 2 r T (m)S 2 = r(— m) which in turn leads to the Onsager relations. For 
the two probe configuration, G(m) = G{— m). For the four probe configuration the 
transmission probabilities satisfy Tjj(m) = tr(tytj-) = m). Focusing on the 

current /volt age configuration: l\ = J 2 = — J 4 , eV\ = fi^—fii and eV^ = — ji2 [29J 
the relation between currents and voltages can be expressed as [30J: 




-a 2 i(m) 




(22) 



where the coefficients ay can be found in Eqs. (4.a-4d) of Ref. [30J. The Onsager 
relations can then be expressed as: 

^•(m) = aji(-m). (23) 

If we choose (say) I\ equal to zero, the relation between the applied current and the spin- 
Hall voltage is: J 2 = Vi(ana 22 — ai2C*2i)/ai2- For phase incoherent conductors, we can 
ignore the interference terms that arise while obtaining the transmission probabilities, 
but the Onsager relations Eq. ( 1231) are unaffected. For a general analysis based on Kubo 
formula see Ref. [32] 

This analysis implies the equivalence of two Hall measurements: (i) setting I\ 
equal to zero and detecting V\ generated by an applied I 2 (Fig. [Tk.) and (ii) switching 
magnetization, setting J 2 equal to zero and detecting V 2 (Fig. [TJd). In other words, 
driving a current 7 2 through the system and detecting the spin Hall voltage with a 
ferromagnetic contact is equivalent to driving a spin accumulation into the SO region 
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Figure 2. a generic four-probe setup for the detection of current-induced spins 



via a ferromagnetic contact that leads to a real Hall voltage detected by normal contacts. 
In the next subsection, we shall exploit this symmetry to gain insight to the boundary 
conditions for a R2DEG|2DEG interface. 

3.2. four-probe setup and boundary conditions 

We now use the Onsager relations from the previous subsection to better understand 
the spin boundary value problem. Consider the four-probe setup in Fig. [IJ When the 
ferromagnetic lead is a Hall contact, the vanishing spin transfer derived by Ref. [26] for a 
(infinitely) wide contact seems to imply that there is neither spin accumulation nor spin 
current near the ferromagnetic reservoir and therefore no Hall voltage. On the other 
hand, in the Onsager equivalent measurement, spins are injected from the ferromagnet 
into the normal region. Since in this case the current is perpendicular to the boundary, 
the spin accumulations can be matched [26] and a spin accumulation in the SO region 
exists. However, the diffusion equation (Q implies that a spin accumulation gives rise 
to a voltage drop in the spin orbit region [3U ITS] . Onsager's relations discussed in the 
previous section imply that these two voltages must be the same provided the injected 
currents are the same. Thus the result for an infinite contact that a current-induced 
spin accumulation can not enter Hall contacts [26] appears to be misleading. In the 
following we shall demonstrate that the spin accumulations around the Hall contact 
must be close up to a numerical factor around the Hall contact. 

We now focus on the current- volt age setup in Fig. [lb. In this case the current is 
perpendicular to the boundary, so the spin accumulations are continuous across an ideal 
R2DEG|2DEG interface. Assuming a diffusive ferromagnet magnetized parallel to the 
current direction and ignoring the resistivity of the normal region, we obtain the spin 
current polarized in the magnetization direction entering the R2DEG: 



where L s = \/Dt s is the (Dyakonov-Perel) spin relaxation length in the R2DEG and 




(24) 



A (m) = L^Dvrto. ■ n + L-^D F v F (l - 5D 2 /4). 



(25) 
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Here, L s p, Dp , vp are the spin relaxation length, diffusion constant and average density 
of states in the ferromagnet, respectively, 5D = {v + D + — v_DJ)/ (upDp), v± and D± are 
the density of states and diffusion constants of the majority and minority spin electrons, 
fj, is a linear function of m of order unity that depends on the details of the geometry 
of the contact. The spin accumulation in the SO region localized within a depth of L s 
at the contact aperture, acts as a dipole source for the diffusion equation: 

V 2 n = V P, (26) 

with dipole density P = —4K s ^ c (z x s)/D. We then estimate the potential drop in the 
Hall direction to be: 

' f drs(r), (27) 



D W 

which is proportional to the integrated spin accumulation 

drs(v) « L S I^-. (28) 



The potential drop is therefore: 

ar I 5D a £ I 5D 



(29) 



° L S W A v F yiTF W A 
up to a numerical constant. 

We now focus on the potential drop in the Onsager-equivalent setting in Fig. UK 
According to the boundary condition Eq. ( [151) . the current induced spin accumulation 
does not enter the normal region. Then the potential drop at the ferromagnet |2DEG 
interface would be zero in contradiction to Onsager's relations. Let us assume that 
the spin accumulations at the R2DEG and 2DEG near the contact are equal to each 
other up to a numerical constant Z, i.e. S2deg = Z sr2Deg- Then the calculation of the 
potential drop proceeds similar to Ref. [32]. Again ignoring the resistance of the 2DEG 
region, we obtain a potential drop as: 

0. = Z- J— -A (30) 

up to a numerical factor. Comparing with Eq. (129]) and noting that we have ignored 
all numerical factors in the calculations above, we conclude that Z must be a numerical 
factor of the order unity in order to satisfy Onsager's relations. In the next section we 
shall consider a model for a narrow contact and show that this is indeed the case. 



4. Model for spin accumulation near a contact 

In this section we focus on the current density and spin accumulation near a finite 
contact between a half-infinite R2DEG and a half- infinite 2DEG (Fig. [3^). The model 
we adopt is sketched in Fig. [3b. Asymptotically, we have a constant current density 
in the left region (R2DEG) in the y direction whereas in the right region (2DEG) the 
charge current density vanishes. The two regions are divided by an infinitely thin and 
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Figure 3. Geometry of the contact: (a) 2D electron gas with a constriction in the 
middle. On the left side there is an applied homogeneous current density which is 
modified near the opening. On the right side, the current density far away from the 
contact as well as the net charge current flowing from the left region to the right region 
is zero. However, there are a finite spin current and a finite spin accumulation in the 
right region. The respective mobilities of the left and right regions are assumed to be 
the same but the Rashba coefficients are different, (b) an idealized version of (a) used 
in the calculations of this section. The origin is chosen at the center of the opening 
with width Wh- 



high potential barrier, except for an opening (the contact) of size Wh centered at (0, 0). 



We note that the solution to this problem closely follows that of an analogous one 
in magnetostatics [35]. We proceed by expressing the chemical potential n in terms of 
the (yet undetermined) solution (ft of the Laplace equation: 

Joy 

vD 

-(f>(x,y) 



n 
n 



+ (f)(x, y) if x < 
if x > 0, 



(31) 



where Jo is the bulk current density in the R2DEG. The asymmetric behaviour of 
in left and right regions is dictated by the current continuity at x = 0. The boundary 
conditions are: 

JoV 



0(0, y) = 

d<j>{0,y) 
dx 



if \y\ < W H /2 
if \y\ > W H /2. 



Next we expand <fi in terms of the modes of the Laplace equation: 



<f){x,y) 



dkA(k)e~ k ^ sm(ky). 



(32) 



(33) 
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The solution to the diffusion equation with the above boundary conditions then reduces 
to that of a dual integral equation: 

dkA(k)sm(ky) = -^ ii\y\<W H /2 

POO 

/ dkkA(k)sm(ky) = if \y\ > W H /2. (34) 
Jo 

Such integral equations arise commonly in potential theory for mixed boundary 
conditions (see Ref. [35J for the solution in 3D). In our case the solution is 

= _*W) (35) 

AvD k 

We may now express the spin accumulations in terms of A(k). For the sake of simplicity, 
we at first disregard the precession term, proportional to K p , in the spin diffusion 
equations Eqs. fl9lfTTl). We shall be particularly interested in the question whether 
current-induced spin accumulation in the spin-orbit coupled region can leak out of 
the contact, into the normal (i.e. no spin-orbit interaction) region. In the bulk of the 
R2DEG, the current is in the y direction, so the current-induced spin accumulation is 
polarized in the x direction. Then the general solution to the spin diffusion equations 
in the R2DEG region is given by: 

Sx(x,y)~ = ( "FT + J + 6s x (x,y), (36) 



uD dy 

where 5s x satisfies the source-free (i.e. zero charge current) diffusion equation that can 
be expanded as: 

POO 

5s x (x,y)= dkB(k)e- Klxl cos(ky), (37) 
Jo 

where k = \Jk 2 + L~ 2 . For the 2DEG side (x > 0), a similar expansion gives: 

POO 

s+(x,y)= dkD(k)e- klxl cos(ky). (38) 



Using the boundary conditions that the spin current is continuous and s x is discontinuous 
by an amount equal to (ar/2)dn/dy |26j, we find that the accumulation in the 2DEG 
satisfies: 

D(k) = -^kA(k) - (K/k)B{k), (39) 
and D(k) is determined from A(k), through the following dual integral equations: 

Warq 2 



[ dqD(q)(l+ J-— ) cos(qy) = - [ dq A(q) 
Jo \ \/q 2 + X 2 J Jo 



cos(qy) (40) 



if \y\ < 1, and 

POO 

dqqD(q) cos(qy) = (41) 



if \y\ > 1. Here we have introduced dimensionless variables q = kWn/2, y = 2y/Wn 
and A = Wh/2L s . In the limit A ^> 1 (wide contact), expanding Eq. (f40l to leading 
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order in A -1 we obtain that D(k) vanishes like A -1 , in agreement with Ref. [26J. In the 
opposite limit A <C 1 (narrow contact), we again expand Eq. f|40l) . this time to leading 
order in A. We then identify the resulting integral equation with the y derivative of 
Eq. (|3lD times ar/2. Thus we show that D(k) = -&-kA(k)/2 solves Eq. flUD up to 
order A 2 corrections. Then the spin accumulation in the 2DEG near a narrow contact 
is given by: 

,+m v) ~ aT dn (°>y) - ^ (A 2 ) 

S *^~ 4 dy - 8uD- [ ] 

We see that the spin accumulation in the 2DEG does not vanish even when the mobilities 
of both sides are equal. For comparison, we also calculate the spin accumulation under 
the assumption that there is no jump in the accumulations. We obtain that in this case 
the spin accumulation is twice as large as s+(0, y). The presence of the term proportional 
to K p generates 2-polarized spin currents going into the 2DEG, owing to the precession 
of y polarized spin accumulation as it diffuses out of the R2DEG, but does not change 
the general picture presented above. We conclude that the choice of the boundary 
condition for spin accumulation near a narrow contact is not important qualitatively, 
because either boundary condition produces identical result up to a numerical factor, in 
agreement with the Onsager's relations. 

5. Numerical results 

In this section, we shall provide a numerical demonstration of the results of the previous 
section, i.e. the possibility of extracting spin accumulations to a normal region with small 
contacts. We focus on the discretized version of the hamiltonian ([1]). Discretization with 
lattice spacing a yields the following tight-binding representation of Ho [36J: 

H ° = 2^2 I 5-^ 4 + f/ ) C ™>™ Cn ' m + Yl ( [ ~~ C n,mCn+l,m (43) 

I n yn ii i ii • 



n,m 



+ H.c. 



where n{m) is the a; (^-coordinate of the site (n, m), a = (ma/h)a. The abbreviation 
c n,m = ( c n,m,+) C n,m-) was used, where cjj ma {cn,m,o) creates(annihilates) an electron at 
site (n, m) with spin orientation a with respect to the z-direction. We also define the 
spin precession length Lgo = na/d, which is related to L s by L so = 2ttL s in the dirty 
limit, but remains well-defined for ballistic systems where there is no spin relaxation. In 
this model, instead of dilute localized scatterers, we shall assume Anderson disorder: the 
dimensionless onsite potential U is set to a different random value U G [—U /2;U /2] 
at each lattice site (n, m) of the disordered region, where Uq accounts for the strength 
of the disorder |43j. The parameter Uq is related to the momentum relaxation rate r 
and the electron mean free path I — v-pr by: 



r = 48a 2 -^, ^ = 48a^, (44) 
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Figure 4. Left panel: Geometry used for numerical calculations: A disordered wire 
with spin orbit coupling and width W connected to two clean leads with spin-orbit 
coupling and to a disordered side-pocket of size 14a x 80a without spin-orbit coupling. 
The colourplot shows the nonequilibrium density (n) averaged over 60000 disorder 
configurations for a system with Lso = 35a, W — 68a, Wh = 20a. The rapid 
oscillations are due to the finite number of channels in the wire. Nevertheless, the 
slow varying part satisfies the diffusion equation. Right panel: Electron density (n) 
of the system shown in the left Panel as a function of vertical coordinate y for fixed 
horizantal coordinate x — 34a (black solid line), x — 65a (red circles) and x = 73a 
(green dashed line). 



where eF = (h 2 /2m*a 2 )~ 1 E? and Ep is the Fermi energy. In the rest of this section, we 
choose Uq — 2 and ep = 0.38 in order to ensure that the transport through the system 
is diffusive. With this choice of parameters the mean free path / « 7.4a is smaller than 
any length scale characterizing the system. 

In order to study the spin accumulation extracted to a normal region we focus on 
the setup shown in Fig. HI where a normal region (i.e. a = 0) with a size of 80a x 14a is 
attached to a Rashba spin-orbit coupled wire of infinite length, width W and constant 
finite spin orbit coupling a > via a contact of size Wh- Disorder of strength Uo is 
present inside the normal region and in the spin-orbit region for —50a < y < 50a. We 
shall use the nonequilibrium Green function method [H] to calculate the lesser Green 
function G < (f;f') which is related to spin accumulation according to 

s x (r)=-~iTr[a x G < (r;r)} (45) 

and to the electron density through 

n(r) = -i Tr [G^f)]. (46) 

Here, we focus on the ensemble averaged accumulations (s x ) and (n). The variances are 
also of interest [401 UT| , but we shall not consider them here. 
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Figure 5. The ratio (s^) / (dn B /dx) as a function of a calculated numerically for two 
different geometries with Wh = 20a and W — 37a (black dots), W — 59a (blue squares) 
and estimated as in Eq. (TiTj) (red line), (s^) and (dn B /dx) have been evaluated by 
averaging over 20000 disorder configurations as well as over the area indicated by the 
blue square shown in the bottom panel of Fig. [5J 



We apply a small bias 5V between the chemical potentials of the top and the bottom 
lead and generate a current in y direction. The left panel of Fig. H] shows the electron 
density (n) inside the system when a current is passed from the top to the bottom. Due 
to the disorder in the central region (—50a < y < 50a) the electron density decreases 
from top to bottom. In the right panel of Fig. H] we show the dependence of (n) on y 
for three different values of x. We observe that (n) decreases linearly in the bulk of the 
spin-orbit region (solid line), showing that the system is diffusive. For x = 65a (circles) 
the side contact at x — 68a disturbs the homogeneous current flow. Inside the normal 
region, x = 73a, (n) is approximately constant (dashed line). 

The current driven by 5V, generates a spin accumulation in the bulk of the R2DEG. 
According to Eq. ffTOl . (s®) = (ar/2)((dn B /dy)) in the bulk. Our simulations agree well 
with the diffusive result as shown in Fig. [5] for large enough a. For smaller values of a, 
Lgo becomes comparable to the overall length of the disorder region L = 100a. In this 
regime ballistic processes can no longer be neglected, causing slight deviations from the 
diffusive theory. 

Having demonstrated that our numerical system is diffusive, we now focus on the 
spin accumulation in the normal region. In Fig. [6], we show the spin density (s x ) averaged 
over 50000 impurity configurations inside three distinct systems with Lgo = 25a. We 
note that in agreement with Ref. [26J, when the interface between R2DEG and 2DEG 
is infinite (top left panel), the spin accumulation in the 2DEG is much smaller than the 
bulk spin accumulation. Nevertheless, when the size of the contact is made smaller 
(top right panel), we observe that the spin accumulation inside the normal region 
increases, reaching a comparable value to the spin accumulation in the bulk when the 
size of the opening is comparable to Lgo (bottom Panel). In order to demonstrate 
this further, we evaluate (s®) by averaging the spin accumulation in the bulk over 
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37a 52a 



Figure 6. Top left panel: Spin accumulation (s x ) in a quantum wire of width W = 52a 
with an abrupt drop of the SO coupling strength at x = 37a from the constant a = tt/25 
{Lso — 25a) for x < 37a to zero on the other side. Top right panel: Spin accumulation 
(s x ) for a system as shown in Fig. |4] with W — 37a, Wjj = 80a and Lso = 25a. Bottom 
panel: Same as top right panel with = 20a. In all three panels, (s x ) is obtained 
by averaging over 50000 disorder configurations. 

the blue square shown in Fig. [6] and (s£) by averaging the spin accumulation in the 
normal conducting side-pocket over the white square shown in Fig. [61 In Fig. [7] we plot 
the ratio (s^.)/(s^) as a function of Lso/Wu, for various values of system and contact 
sizes. We observe that starting from small Lso/Wh, the spin accumulation increases 
with L S o/Wh, approaching to ~ 0.5 — 0.7. This value is in between the estimates 
0.5 and 1.0 based on diffusion equations using the boundary conditions Eq. (1151) and 
Eq. ([ED respectively. For small values of L so /Wu (Fig. [3, left panel), (s^)/(s^) is of 
order (Lso/Wu) in agreement with the analytical calculation above. We note, however, 
that in this limit the system we considered is close to the clean limit L$o ~ ^? where 
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Figure 7. Left Panel: Average spin accumulation inside the normal region (s£) 
relative to the accumulation in the bulk of the spin-orbit region (sjj?) for various 
geometries averaged over 20000 disorder configurations as a function of Lso/Wh. 
Right Panel: (s£)/(s^) for two different geometries averaged over 60000 disorder 
configurations. 



deviations from the diffusion equations might be expected. Currently, we are working 
on larger systems in order to explore small Lso/Wh hi the dirty limit [35] . 

6. Conclusions 

In this work, we considered the problem of extracting current-induced spins generated 
in a region with spin-orbit coupling into a region with vanishing (or small) spin-orbit 
coupling, where the spin relaxation time is long. To this end we focused on the spin 
boundary conditions between a spin-orbit coupled region and a normal region. Although 
for an infinite interface the spins are confined to the spin-orbit region via the boundary 
spin Hall effect, we have shown by solving a model problem as well as doing numerical 
simulations that for a finite interface the spin accumulations generated in the spin-orbit 
region can be extracted to a normal region. The amount of extracted spin accumulation 
is equal to that of the spin-orbit region up to a geometrical factor of order unity. 
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